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L, F. Guseman, Jr. and Bruce P. Marlon 


1. Introduction 


This report presents preliminary numerical results obtained from the 
application of a linear feature selection technique to the determination of 
combinations of passes which best discriminate between a given set of 
crops in a given area of interest. The results obtained are not purported 
to hold in a general situation, but only for the given set of crops and 
the given (but unknown) levels of several factors (such as soil type, 
fertilizer practice, etc.) holding in the area of interest. Nevertheless, 
some insight into optimal choices of passes can be gained through the 
application of the technique to data sets with different levels of those 
factors which are known to affect the spectral signatures. The technique 
could also be useful in determining those factors which significantly 
affect the choice of passes used for discrimination between a given set 
of crops. 

For a given set of r registered ERTS passes, the input data consists 
of a Ar-dimensional mean vector and a 4r x 4r covariance matrix for 
each of the m classes (assumed to be multivariate normal). For each of 
the combinations of r passes taken k at a time (1 < k < r) the portions 
of the mean vectors and covariance matrices pertaining to the k passes 
i ,i ,.,.,i. are used to compute a 4k-dlmensional vector B(i ,...,!,) 
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of norm one whose linear combination of the channels in the passes 
i ,i , ...» i. minimizes the probability of misclassif ication, g(B), for 
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the transformed (one-dimensional) densities of the m classes. The 
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2. Numerical Results-- Hill County North 

Input data for the runs described in this section consisted of 
24-dimensional mean vectors and 24 x 24 covariance matrices computed from 
training data using six registered passes of ERTS data from Hill County 
North. The six registered ERTS passes (Identified below) were, for 


purposes of this note, assigned the following pass numbers and channel 


numbers. 

given. 

The associated 

biological phases for Winter Wheat 

are also 

Pass 

Channel 

Date 

Scene 

Biological 

Number 

Numbers 

of Pass 

Number 

Phase 

l (Winter Wheat) 

1 

1-4 

Aug. 20, 1973 

1393-17392 

Harvest 

2 

5-8 

Aug. 3, 1973 

1376-17452 

Yellow/ Harvest 

3 

9-12 

July 16, 1973 

1358-17453 

Heading/Yellow 

4 

13-16 

June 27, 1973 

1339-17400 

Heading 

5 

17-20 

June 11, 1973 

1332-17460 

Heading 

6 

21-24 

May 23, 1973 

1304-17461 

Jointing (Greening) 


(Note that the passes are numbered in Reverse chronological order) 


The crops under consideration and their designated class numbers 

were: 


Class 


Crop 


1 

2 

3, 

4 

5 

6 


Barley 

Stubble 

Grass 

Winter Wheat 
Spring Wheat 
Fallow 


For each of the sixty-three combinations of 6 passes taken k at a 
time (1 < k < 6), the appropriate signatures were used to compute a 
vector B of norm one (and appropriate dimension) whose linear combination 
of the channels minimizes the one-dimensional probability of misclassifi- 
cation. The resulting values, g(B), of the one-dimensional probability 
of mi8classif ication determined by the optimal B vectors for the various 
combinations of passes appear in Tables 1-5 below. The value for g(B) 
using all channels within all bIx passes was .2699. The graphs of the 
six transformed density functions determined by the optimal B vector 
for the optimal one, two, three, four, five and six passes appear in 
Figures 1-6, respectively. 

Additional runs were made using all six passes wh'Ve Spring Wheat 
and Winter Wheat were alternately deleted and the remaining five crops 
(renumbered in same order) considered. The results from these runs 
appear in Figures 7 and 8, respectively. 

The confusion matrices associated with Figures 1 $ -appear in 


Tables 6 and 7. 

























1.2, 3,4 

1.2, 3, 5 

1.2, 3, 6 

1.2 4,5 

1.2, 4, 6 

g(B) 

.3917 

.3213 

.3226 

a 

.3247 

Passes 

1.2, 5, 6 

1.3, 4, 5 

1. 3.4,6 

SB 

1.4, 5,6 

t(B) 

.2764 

.3377 

.3521 

.2873 

.2873 

Passes 

2, 3, 4, 5 

2, 3, 4, 6 

2, 3, 5, 6 

2,4, 5,6 

3, 4, 5, 6 

g(B) 

.3229 

.3738 

.2925 

.2909 



Table 4. 


Values of g(B) 



Passes 

1.2, 3, 4. 5 

1,2, 3, 4, 6 

1,2, 3, 5, 6 

1,2, 4, 5, 6 

1.3, 4, 5, 6 

2, 3,4. 5,6 

g (B) 

.3171 

.3140 

.2733 

.2718 

.2842 

.2886 


Table 5. 


Values of g(B) 




Optr^sl B, Transformed densities: Best one pass- 








Fig. 3. Optimal B, g(B), Transformed densities: Best three passes-1,5.6. 





Channels 








Channels 1-8, 13-24 
B( 1) - -.0244 
B( 2) - -.0678 

H.BHE-I2 A ~ ?! ?? " ■*???? 






Optimal B, g(B), Transformed densities: All passes-1 ,2, 3, 4, 5, 6 



lasses (no Winter Wheat): Optimal B„ g(B), Transformed densities: All passes 










Fig. 8. Five classes (no Spring Wheat): Optimal B, g(B), Transformed densities: All passes. 


668 .000 .000 .003 .050 .000 .525 .000 .000 .005 .121 



Table 6. Confusion matrices 


493 .000 .047 .003 .168 .000 / .528 .000 .047 .0C2 .176 .001 
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3 . Rema rks 

The preliminary numerical results in the previous section Indicate 
that fewer optimally chosen passes could be used while maintaining most 
ol the dlscrlmatory power Inherently available in the original data. 

Again, we note that the results give Insight only for the crops con- 
sidered -?nd levels of factors reflected by the data used. However, by 
Identifying the various factors (an levels thereof) affecting the 
spectral signatures, and by formulating a regression model one could 
use the feature selection technique to determine the regression 
coefficients for predicting optimal passes for a given set of crops. 

Another use of the feature selection technique as applied to 
multiple-pass registered data is the generation of enhanced grey 
scale displays by using a single linear combination of all channels of all 
designated passes as opposed to a single channel within a single pa'-.s. 




